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Quantifier 



Data members 

(required) 

Array dimensions J, J, K 
Scalar type 

Density array (address) 
Method-choice flags 

(as needed) 
Bound_density array 
Free_density array 
Binding_rate array 
Creationjrate array 
IxJxK segmentation array (address) 
Segment-surface list (address) 
Transport model (address) 
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Function members 

(scalar-valued) 

Volume_in_range( [Co], [G], [segment] ) 
lntegrate( [segment | triangulated_surface] ) 
Efflux( triangulated_surface ) 

(array-valued) 

Difference( subtrahend, [segment], [divisor], [target] ) 
MuItMmage( arrayjist, algorithm, [target] ) 
lmage_process( operator, [segment], [target] ) 



+ member-setting functions 

J 



64<T 
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TemporalQuantifier 



711 ] Data members 
(required) 

s Array dimensions /, J, K 
s NumberiVof times 
^Times /i, f2, ... , <w 
Scalar type 
Density array list 
Method-choice flags 

(as needed) 

Sound density array list 
-Free density array list 
-Binding rate array list 
-Creation rate array list 
JxJxK segmentation array (static or list) 
-Segment-surface (static or list) 
-Transport model (address) 
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Function members 

(scalar-valued) 
s Volume_in_range( [Co], [G], [segment], [times]) 
lntegrate( [segment | triangulated_surface], [times] ) 
"Efflux( triangulated_surface_list, [times] ) 

(array-valued) 

Rate_of_change( [segment], [divisor], [target], [times] ) 
-Multi_image( arrayjist, algorithm, [target], [times] ) 
Jmage_process( operator, [segment], [target], [times] ) 
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+ member-setting functions 



J. 
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Initialize pointer on 
first triangle. 
Volume6 = 0.0 



810 



Volume += prism volume 



last triangle? 



820/-\_y 
82 K / 



No 



next 
triangle 



Yes 
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Return volume6/ 
6.0 
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911 
912 

920 
921 
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Initialize pointer on 
first gridpoint. 
Count = 0 . 



Ci< c(point)< C 2 ? 



Yes 



In segment S? 



"V 
Yes 
T 



No 



Count += 1 



No 



v y 



Last grid point? 



No 



next 
point 



Yes 



Return count xH 
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1050 



Initialize pointer on first 
volume element. 

V = 0.0 



Find 8-corner type T with respect to 
above/below threshold C 



Interpolate edge lengths relevant to T 



V += polynomial relevant to T 




Last element? 



No 



next 
element 



Yes 




Ficofte ID 



FV 
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1030 



1140 
1141 



Initialize 'perhaps in S' 
volume element iterator. 

V = 0.0 



FV 



All corners outside segment S? 



All corners Inside S? 




Yes 



Yes 



V += fractional content 



Find 8-corner type T with respect to 
condition 'above/below threshold C 



Interpolate edge lengths relevant to T 



V += polynomial relevant to T 
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1120 



,1130 



Last element? 



No 

T 



next 
element 



Yes 



Return VxH 
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Input array parameters, 

range limits &, C 2 , 
optionally segment S. 



Set range 
condition 

<C2 



FV or FV 



1 

Set r 
cond 
< 


ange 
ition 

C 2 






FV or FV 




Return difference 
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Estimate dc/dt over grid, 
in array with address D 



Select a boundary B where 
c is negligible 



Fix array for free-density values, 
with address A. 



Pass D, A and B to transport 
model inversion routine. 



Return A 
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Create array e, 
initialize values to 0.0 



Triangle index i = 0 
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Compute normal vector n of triangle i 



Compute mid-point w of triangle i 



Pass w to transport model 



Receive vector F from transport model 



e[N] += n. F, e[i] = n • F/(triangle area) 
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Last triangle? 



No 



index i += 1 



Yes 



Return A 
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1430 



1440 



1441 



1450 
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